Stellar weak decay rates in neutron-deficient medium-mass nuclei 
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Weak decay rates under stellar density and temperature conditions holding at the rapid proton 
capture process are studied in neutron-deficient medium-mass waiting point nuclei extending from 
Ni up to Sn. Neighboring isotopes to these waiting point nuclei are also included in the analysis. The 
nuclear structure part of the problem is described within a deformed Skyrme Hartree-Fock -I- BCS -f 
QRPA approach, which reproduces not only the beta-decay half-lives but also the available Gamow- 
Teller strength distributions, measured under terrestrial conditions. The various sensitivities of the 
decay rates to both density and temperature are discussed. In particular, we study the impact of 
contributions coming from thermally populated excited states in the parent nucleus, as well as the 
competition between beta decays and continuum electron captures. 

PACS numbers: 23.40.-s,21.60. Jz,26.30.Ca,27.50.+e 



I. INTRODUCTION 

An accurate understanding of most astrophysical pro- 
cesses requires necessarily information from nuclear 
physics, which provides the input to deal with network 
calculations and astrophysical simulations (see [l|, and 
references therein) . Obviously, nuclear physics uncertain- 
ties will finally affect the reliability of the description of 
those astrophysical processes. This is especially relevant 
in the case of explosive phenomena, which involve knowl- 
edge of the properties of exotic nuclei, not well explored 
yet. Thus, most of the astrophysical simulations of these 
violent events must be built on nuclear-model predictions 
of limited quality and accuracy. This is in particular the 
case of the X-ray bursts (XRBs) 0-1^, which are gen- 
erated by a thermonuclear runaway in the hydrogen-rich 
environment of an accreting neutron star that is fed from 
a red giant binary companion close enough to allow for 
mass transfer. 

Type I XRBs are typically characterized by a rapid 
increase in luminosity generating burst energies of 10'^^ — 
10"'" ergs, which are typically a factor 100 larger than 
the steady luminosity. The luminosity suffers a sharp 
raise of about 1 — 10 s followed by a gradual softening 
with time scales between 10 and 100 s. These bursts are 
recurrent with time scales ranging from hours to days. 
The properties of XRBs are particularly dependent on 
the accretion rate. Typical accretion rates for type I 
XRBs are about 10~* — 10~^Mq yr~^. Lower accretion 
rates lead to weaker flashes while larger accretion rates 
lead to stable burning on the surface of the neutron star. 

The ignition of XRBs takes place when the tempera- 
ture (T) and the density [p) in the accreted disk become 
high enough to allow a breakout from the hot CNO cycle. 
Peak conditions of T = 1 - 3 GK and p = lO*' - lO"^ g 
cm^'^ are reached and eventually, this scenario allows the 
development of the nucleosynthesis rapid proton capture 
{rp) process which is characterized by proton cap- 

ture reaction rates that are orders of magnitude faster 
than any other competing process, in particular /3-decay. 
It produces rapid nucleosynthesis on the proton-rich side 



of stability toward heavier proton-rich nuclei reaching nu- 
clei with A > 100, as it have been studied in Ref. [3], 
where it was shown that the rp process ends in a closed 
SnSbTe cycle. It also explains the energy and luminosity 
profiles observed in XRBs. 

Nuclear reaction network calculations, which may in- 
volve as much as several thousand nuclear processes, are 
performed to follow the time evolution of the isotopic 
abundances, to determine the amount of energy released 
by nuclear reactions, and to find the reaction path for 
the rp process jsl-fioj. In general, the reaction path fol- 
lows a series of fast proton-capture reactions until the 
dripline is reached and further proton capture is inhib- 
ited by a strong reverse photodisintegration reaction. At 
this point, the process may only proceed through a beta 
decay or a less probable double proton capture. Then the 
reaction flow has to wait for a relatively slow ^S-decay and 
the respective nucleus is called a waiting point (WP). The 
short time scale of the rp process (around 100 s) makes 
highly significant any mechanism that may affect the pro- 
cess in some seconds and the half-lives of the WP nuclei 
are of this order. Therefore, the half-lives of the WP 
nuclei along the reaction path determine the time scale 
of the nucleosynthesis process and the produced isotopic 
abundances. In this respect, the weak decay rates of 
neutron-deficient medium-mass nuclei under stellar con- 
ditions play a relevant role to understand the rp process. 

Although the products of the nucleosynthesis rp pro- 
cess are not expected to be ejected from type I XRBs 
due to the strength of the neutron star gravitational field, 
there are other speculative sites for the occurrence of rp 
processes. This is the case of core collapse supernovae 
that might supply suitable physical conditions for the rp 
process provided neutrino- induced reactions are included 
in the nucleosynthesis calculations [ll|. These reactions 
have to be included to bypass the slow beta decays at the 
WP nuclei via capture reactions of neutrons, which are 
created from the antielectron neutrino absorption by free 
protons [T2|. Contrary to the XRBs, these scenarios will 
finally lead to the ejection of the nucleosynthetic products 
and thus contribute to the galactic chemical evolution. 

Since the pioneering work of Fuller, Fowler and New- 
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man , where the general formahsm to calculate weak- 
interaction rates in stellar environments as a function of 
density and temperature was introduced, improvements 
have been focused on the description of the nuclear struc- 
ture aspect of the problem. Different approaches to de- 
scribe the nuclear structure involved in the stellar weak 
decay rates can be found in the literature. They are ba- 
sically divided into Shell Model [l^, [1 3 or quasiparticle 
random phase approximation (QRPA) [l5l - [l7| categories. 
Certainly, the nuclear structure problem involved in the 
calculation of these rates must be treated in a reliable 
way. In particular, this implies that the nuclear mod- 
els should be able to describe at least the experimental 
information available on the decay properties (Gamow- 
Teller strength distributions and /3-decay half-lives) mea- 
sured under terrestrial conditions. Although these decay 
properties may be different at the high p and T existing 
in rp process scenarios, success in describing the decay 
properties in terrestrial conditions is a requirement for a 
reliable calculation of the weak decay rates in more gen- 
eral conditions. With this aim in mind, we study here 
the dependence of the decay rates on both p and T using 
a QRPA approach based on a selfconsistent deformed 
Hartree-Fock (HF) mean field. Deformation has to be 
taken into account because the reaction path in the rp 
process crosses a region of highly deformed nuclei around 
A — 70 — 80. This nuclear model has been tested suc- 
cessfully (see [l^ and references therein) and reproduces 
very reasonably the experimental information available 
on both bulk and decay properties of medium-mass nu- 
clei. In this work we focus our attention to the even-even 
WP Ni, Zn, Ge, Se, Kr, Sr, Zr, Mo, Ru, Pd, Cd, and Sn 
isotopes, as well as to their closer even-even neighbors. 

The paper is organized as follows. In Section |lT] the 
weak decay rates are introduced as functions of density 
and temperature and their nuclear structure and phase 
space components are studied. Section IIIII contains the 
results. First, we study the decay properties under terres- 
trial conditions, and secondly as functions of both den- 
sities and temperatures at the rp process. Section IIVI 
contains the conclusions of this work. 



II. WEAK DECAY RATES 

There are several distinctions between terrestrial and 
stellar decay rates caused by the effect of high p and 
T. The main effect of T is directly related to the ther- 
mal population of excited states in the decaying nucleus, 
accompanied by the corresponding depopulation of the 
ground states. The weak-decay rates of excited states can 
be significantly different from those of the ground state 
and a case by case consideration is needed. Another ef- 
fect related to the high p and T comes from the fact that 
atoms in these scenarios are completely ionized and con- 
sequently electrons are no longer bound to the nuclei, but 
forming a degenerate plasma obeying a Fermi-Dirac dis- 
tribution. This opens the possibility for continuum elec- 



tron capture {cEC), in contrast to the orbital electron 
capture (oEC) produced by bound electrons in the atom 
under terrestrial conditions. These effects make weak in- 
teraction rates in the stellar interior sensitive functions 
of T and p, with T = 1.5 GK and p = 10^ g cm'^, as the 
most significant conditions for the rp process 
The decay rate of the parent nucleus is given by 

^^Y^^^2jh+l^-E.n,T)^ (1) 

i 

where G ^ J2i i'^Ji + 1) e--^'/^^^) is the partition func- 
tion, Ji{Ei) is the angular momentum (excitation energy) 
of the parent nucleus state i, and thermal equilibrium is 
assumed. In principle, the sum extends over all popu- 
lated states in the parent nucleus up to the proton separa- 
tion energy. However, since the range of temperatures for 
the rp process peaks at T = 1.5 GK {kT ^ 300 keV), only 
a few low-lying excited states are expected to contribute 
significantly in the decay. Specifically, we consider in 
this work all the collective low-lying excited states be- 
low 1 MeV [l^. Two-quasiparticle excitations in even- 
even nuclei will appear at an excitation energy above 2 
MeV, which is a typical energy to break a pair in these 
isotopes. Hence, they can be safely neglected at these 
temperatures. As an example, the maximum population 
appears for the lowest of these states {E2+ =261 keV in 
■^^Sr), which at T=1.5 GK is 12%, while the ground state 
still contributes with 88%. 

The decay rate for the parent state i is given by 

/ 

where the sum extends over all the states in the final 
nucleus reached in the decay process. The rate A;/ from 
the initial state i to the final state / is given by 

A./ = ^B,/$./(p,T), (3) 

where D — 6146 s. This expression is decomposed into a 
nuclear structure part Bif that contains the transition 
probabilities for allowed Fermi (F) and Gamow- Teller 
(GT) transitions, 

B,f^B,f{GT)+B,fiF), (4) 

and a phase space factor ^if, which is a sensitive function 
of p and T. The theoretical description of both Bif and 
^if are explained in the next subsections. 

A. Nuclear Structure 

The nuclear structure part of the problem is described 
within the QRPA formalism. Various approaches have 
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been developed in the past to describe the spin-isospin 
nuclear excitations in QRPA (20l - [3H . In this subsection 
we show briefly the theoretical framework used in this pa- 
per to describe the nuclear part of the decay rates in the 
neutron-deficient nuclei considered in this work. More 
details of the formalism can be found in Refs. [29l - l3TI |. 

The method starts with a self-consistent deformed 
Hartree-Fock mean field formalism obtained with Skyrme 
interactions, including pairing correlations. The single- 
particle energies, wave functions, and occupation proba- 
bilities are generated from this mean field. In this work 
we have chosen the Skyrme force SLy4 [s^ as a repre- 
sentative of the Skyrme forces. This particular force in- 
cludes some selected properties of unstable nuclei in the 
adjusting procedure of the parameters. It is one of the 
most successful Skyrme forces and has been extensively 
studied in the last years. 

The solution of the HF equation is found by using the 
formalism developed in Ref. 33] , assuming time reversal 
and axial symmetry. The single-particle wave functions 
are expanded in terms of the eigenstates of an axially 
symmetric harmonic oscillator in cylindrical coordinates, 
using twelve major shells. The method also includes pair- 
ing between like nucleons in BCS approximation with 
fixed gap parameters for protons and neutrons, which are 
determined phenomenologically from the odd-even mass 
differences involving the experimental binding energies 

M- 

The potential energy curves are analyzed as a function 
of the quadrupole deformation /3, 



Qc 



5 A{r^ 



(5) 



written in terms of the mass quadrupole moment Qo and 
the mean square radius (r^). For that purpose, con- 
strained HF calculations are performed with a quadratic 
constraint (ssj . The HF energy is minimized under the 
constraint of keeping fixed the nuclear deformation. Cal- 
culations for GT strengths are performed subsequently 
for the various equilibrium shapes of each nucleus, that 
is, for the solutions, in general deformed, for which min- 
ima are obtained in the energy curves. Since decays con- 
necting different shapes are disfavored, similar shapes are 
assumed for the ground state of the parent nucleus and 
for all populated states in the daughter nucleus. The 
validity of this assumption was discussed for example in 
Refs. Illll- 

To describe GT transitions, a spin-isospin residual in- 
teraction is added to the mean field and treated in a de- 
formed proton-neutron QRPA. This interaction contains 
two parts, a particle-hole {ph) and a particle-particle 
(pp) . The interaction in the ph channel is responsible for 
the position and structure of the GT resonance [l^, 
and it can be derived consistently from the same Skyrme 
interaction used to generate the mean field, through the 
second derivatives of the energy density functional with 
respect to the one-body densities. The ph residual inter- 



action is finally expressed in a separable form by averag- 
ing the resulting contact force over the nuclear volume 
[29| . The pp part is a neutron-proton pairing force in the 
= coupling channel, which is also introduced as 
a separable force i23|, i30|] . The strength of the pp resid- 
ual interaction in this theoretical approach is not derived 
self-consistently from the SLy4 force used to obtain the 
mean field, but nevertheless it has been fixed in accor- 
dance to it. This strength is usually fitted to reproduce 
globally the experimental half-lives. Various attempts 
have been done in the past to fix this strength 2J], ar- 
riving to expressions that depend on the model used to 
describe the mean field, Nilsson model in the above refer- 
ence. In previous works [13, H^-fUl we have studied the 
sensitivity of the GT strength distributions to the various 
ingredients contributing to the deformed QRPA-like cal- 
culations, namely to the nucleon-nucleon effective force, 
to pairing correlations, and to residual interactions. We 
found different sensitivities to them. In this work, all 
of these ingredients have been fixed to the most reason- 
able choices found previously [11] and mentioned above. 
In particular we use the coupling strengths Xgt =0.15 
MeV and k^Jj, = 0.03 MeV. 

The proton-neutron QRPA phonon operator for GT 
excitations in even-even nuclei is written as 



r+ 



E 



(6) 



where a"*" (a) are quasiparticle creation (annihilation) 
operators, luk are the QRPA excitation energies, and 
X^^ ,Y^J'^ the forward and backward amplitudes, re- 
spectively. For even-even nuclei the allowed GT tran- 
sition amplitudes in the intrinsic frame connecting the 
QRPA ground state |0) (r^^^ |0) = 0) to one-phonon 
states \ujk} {T^,, |0) = \ujk}), are given by 



(7) 



where 



with 



M-" = 5](g..x;j'-H-g,,y--), (8) 

MX" = , (9) 



77riy = u^v-^T."!^, q^r^ = v^u^'E'^i^, (10) 



v's are occupation amplitudes {u^ = 1 — w^) and spin 
matrix elements connecting neutron and proton states 
with spin operators 



(11) 



The GT strength for a transition from an initial state 
i to a final state / is given by 
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2J^ 



yy / off j 



(12) 



where {gA/9v)cff = 0.74(g^/gv)baro is an effective 
quenched value. For the transition IiKi{0^0) — >• 
IfKf{l^K) in the laboratory system, the energy dis- 
tribution of the GT strength Buj{GT^) is expressed in 
terms of the intrinsic amplitudes in Eq. ([7]) as 



gA 
ffV/cff 



Kfi 



+2{u:K\ait^\Qf 5ka 



(13) 



To obtain this expression, the initial and final states in 
the laboratory frame have been expressed in terms of 
the intrinsic states using the Bohr-Mottelson factoriza- 
tion [40|]. 

Concerning Fermi transitions, the Fermi operator is 
the isospin ladder operator T±, which commutes with 
the nuclear part of the Hamiltonian excluding the small 
Coulomb component. Then, superallowed Fermi tran- 
sitions (0^ — O"'') only occur between members of an 
isospin multiplet. The Fermi strength is narrowly con- 
centrated in the isobaric analog state (IAS) of the ground 
state of the decaying nucleus. Thus, neglecting effects 
from isospin mixing one has 



B. Phase Space Factors 

The phase space factor contains two components, elec- 
tron capture [EG) and /3+ decay 



= + . (17) 

In the case of (3~^ /EG decay in the laboratory, EG 
arises from orbital electrons in the atom and the phase 
space factor is given by [i^l 



ff,oEC _ „2^2 n 



(18) 



where x denotes the atomic subshell from which the elec- 
tron is captured, q is the neutrino energy, g is the radial 
component of the bound-state electron wave function at 
the nucleus, and B stands for other exchange and overlap 
corrections [43 |. 

In rp-process stellar scenarios, the phase space factor 
for cEG is given by 



^cEC 



ujp{Q,f+ uj)^F{Z, oj) 
xS,-{uj)[l-S,{Q,f+oj)]dij. (19) 



The phase space factor for positron emission pro- 
cess is given by 



1 



2.h + l 



3 



if ^T{T+l)-T,J, 



(14) 

where T is the nuclear isospin and Tz — {N — Z)/2 its 
third component. The Bif{F^) strength of our concern 
here reduces to B{F+) = {Z - N) = 2 for the iT,Tz) ^ 
(1, -1) isotopes in the decay (Z,A^) (Z-1, -1-1) with 
Z = N + 2. For these transitions the excitation energy 
of the IAS in the daughter nucleus is given by [1, [lOj 

Bias = (ME), ~ (ME) f + 0.7824 - AEc MeV, (15) 

where ME is the atomic mass excess. The Coulomb dis- 
placement energy AEc between pairs of isobaric analog 
levels is given by 



AEc = 1.4144Z/A 



1/3 



0.9127MeV, 



(16) 



where Z = {Zi + Zf)/2. This expression was obtained in 



Ref. [411 from a fitting to data corresponding to levels 
with isospin T = 1. In any case, Fermi transitions are 
only important for the /3+ decay of neutron-deficient light 
nuclei with Z > N {Tz < 0), where the IAS can be 
reached energetically. Thus, although they have been 
considered in the calculations of the terrestrial half-lives, 
only the dominant GT transitions are included in the 
stellar decay rates. 



if ~ 



ujp{Q^f-ufF{^Z+\,Lo) 



x[l-S,^{uj)][l-S,{Q,f-Lo)]duj. (20) 



In these expressions lo is the total energy of the positron 
in TTieC^ units, p = — 1 is the momentum in TOeC 
units, and Qif is the total energy available in m^c^ units 



{Mp- Md + E,- Ef) 



(21) 



which is is written in terms of the nuclear masses of par- 
ent (Mp) and daughter (M^;) nuclei and their excitation 
energies Ei and Ef, respectively. F{Z,uj) is the Fermi 
function [l^ that takes into account the distortion of the 
/3-particle wave function due to the Coulomb interaction. 



F{Z,uj) = 2(1 + 7)(2pi?)-^(i-^)e"^ iwl^,'!!!' ' (22) 



[r(27 + l)] = 



where 7 — -^/l — (aZ)^ ; y ~ aZcu/p ; a is the fine 
structure constant and R the nuclear radius. The lower 
integration limit in the cEG expression is given by = 1 
if Qif > -1, or W£ = \Qif \ if Qif < -1. 

S^-, Se+, and S'y, are the electron, positron, and neu- 
trino distribution functions, respectively. Its presence 
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inhibits or enhances the phase space available. In rp 
scenarios the commonly accepted assumptions @ state 
that — since neutrinos and antineutrinos can es- 
cape freely from the interior of the star and then they do 
not block the emission of these particles in the capture 
or decay processes. Positron distributions become only 
important at higher T {kT > 1 MeV) when positrons 
appear via pair creation, but at the temperatures con- 
sidered here we take 3^+ = 0. The electron distribution 
is described as a Fermi-Dirac distribution 



III. RESULTS FOR WEAK DECAY RATES 

In this section we present first the results for the po- 
tential energy curves. Then, we show the results for the 
decay properties, GT strength distributions and /3-decay 
half-lives, under terrestrial conditions comparing them 
with the available experimental information. Finally, we 
present the results for the stellar weak decay rates un- 
der density and temperature conditions implied in the rp 
process. 



' exp [(^^ - /ie) /(fcT)] + 1 ' ^ ' 

assuming that nuclei at these temperatures are fully ion- 
ized and electrons are not bound to nuclei. The chemical 
potential /ig is determined from the expression 

py^-^,{^yfiS.-S^.)pdp, (24) 

in (mol/cm'^) units, p is the baryon density (g/cm'^), 
is the electron-to-baryon ratio (mol/g), and Na is Avo- 
gadro's number (mol~^). 

Under the assumptions 3^+ = = mentioned 
above, the phase space factors for /S"*" decay in Eq. (PO)) 
are independent of the density and temperature. The 
only dependence of the decay rates on T arises from 
the thermal population of excited parent states. On the 
other hand, the phase space factor for cEC in Eq. 
is a function of both pY^ and T, through the electron 
distribution S^- ■ 

The phase space factors increase with Qi f and thus the 
decay rates are more sensitive to the strength Bi f located 
at low excitation energies of the daughter nucleus. It is 
also interesting to notice the relative importance of both 
decay and electron capture phase space factors (see 
Fig. 3 in Ref . [iTj ) . In general, the former dominates 
at sufficiently high Qif (low excitation energies in the 
daughter nucleus), while the latter is always dominant 
at sufficiently low Qif (high excitation energies in the 
daughter nucleus). 

The /3-decay half-life in the laboratory is obtained by 
summing all the allowed transition strengths to states in 
the daughter nucleus with excitation energies lying below 
the corresponding Qec energy, and weighted with the 
phase space factors. 



Tul = ]^2 = i5 S [S./(GT) + B.,(F)]$f;/-^ 

0<Ej<Qec 

(25) 

where the Qec energy is given by 



Qec = Mp - Md + me^ Qfj+ + 2me . (26) 



A. Potential Energy Curves 

In Fig. [1] we can see the potential energy curves for 
the even-even Ni, Zn, Ge, Se, Kr, Sr, Zr, Mo, Ru, Pd, 
Cd, and Sn nuclei in the vicinity of the N — Z isotopes 
considered in this work. We show the energies relative 
to that of the ground state plotted as a function of the 
quadrupole deformation j3 in Eq. ([S]). They are obtained 
from constrained HF-I-BCS calculations with the Skyrme 
force SLy4 

The nuclei studied here cover a whole proton shell 
ranging from magic number Z = 28 (Ni isotopes) up 
to magic number Z = 50 (Sn isotopes). The isotopes 
considered are the predicted WP nuclei, which in most 
cases correspond to N = Z, and their neighbor isotopes. 

Then, it is expected that the lighter and heavier nu- 
clei close to Z = 28 and Z — 50, respectively, have a 
tendency to be spherical. The spherical shapes in these 
isotopes show sharply peaked profile that become shallow 
minima as one moves away from Z — 28 or Z = 50, and 
finally deformed shapes are developed as one approaches 
mid-shell nuclei. The profiles of the latter exhibit a rich 
structure giving raise to shape coexistence when various 
minima at close energies are located at different defor- 
mations. 

It is also worth mentioning the correlations observed 
between mirror nuclei interchanging the number of neu- 
trons and protons. Thus, we see the remarkable similar- 
ity between the profiles of ^^Ge (Z = 32,N = 34) and 
66Se {Z = 34,7V = 32), between ^^Se {Z = 34, = 36) 
and ™Kr {Z = 36,N ^ 34), and between ^^Ki [Z = 
36, N = 38) and ^^Sr {Z = 38, TV = 36). 

These results are in qualitative agreement with similar 
ones obtained in this mass region from different theoreti- 
cal approaches. Just to give some examples, shape tran- 
sition and shape coexistence were discussed in A ~ 80 
nuclei within a configuration-dependent shell-correction 
approach based on a deformed Woods-Saxon potential 
[43|. Relativistic mean field calculations in this mass re- 
gion have also been reported in Ref. [3| • Nonrelativistic 
calculations are also available from both Skyrme [45|f47| 
and Gogny [4^ forces, as well as from the complex VAM- 
PIR approach 

Experimental evidence of shape coexistence in this 
mass region has become available in the last years 
l65j . and by now this is a well established characteristic 
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feature in the neutron-deficient ^ = 70 — 80 mass region. 



B. Laboratory Gamow- Teller strength and 
half-lives 

While the half-lives give only a limited information of 
the decay (different strength distributions may lead to 
the same half-life), the strength distribution contains all 
the information. It is of great interest to study the decay 
rates under stellar rp conditions using a nuclear struc- 
ture model that reproduces the strength distributions 
and half-lives under terrestrial conditions. 

In the next figures, we show the results obtained for the 
energy distributions of the GT strength corresponding to 
the equilibrium shapes for which we obtained minima in 
the potential energy curves in Fig. [TJ The GT strength 
is plotted versus the excitation energy of the daughter 
nucleus E^^^ = Ef (MeV). 

Fig. [2] ([3]) contains the results for the isotopes Ni, 
Zn, Ge, Se, Kr, and Sr (Zr, Mo, Ru, Pd, Cd, and Sn). 
We show the energy distributions of the individual GT 
strengths in the case of the ground state shapes. We 
also show the continuous distributions for both ground 
state and possible shape isomers, obtained by folding the 
strength with 1 MeV width Breit-Wigner functions. The 
vertical arrows show the Q ec energy, as well as the pro- 
ton separation energy in the daughter nucleus, both taken 
from experiment [34| . 

It is worth noticing that in general both deformations 
produce quite similar GT strength distributions on a 
global scale. The main exceptions correspond to the com- 
parison between spherical and deformed shapes, where 
clear differences can be observed. In any case, the small 
differences among the various shapes at the low energy 
tails (below the Qec) of the GT strength distributions 
lead to sizable effects in the /3-decay half-lives. These 
differences can be better seen because of the logarithmic 
scale. 

Experimental information on GT strength distribu- 
tions are mainly available for ^^Kr ""^Kr [g^], ''^Sr 
[il], and ^°^^^"^Sn [69] isotopes, where /J^-decay exper- 
iments have been performed with total absorption spec- 
troscopy techniques, allowing the extraction of the GT 
strength in practically the whole Q-energy window. In 
Ref. [l^ a comparison between similar calculations to 
those in this work and the experimental data for Kr and 
Sr isotopes was carried out. In general, good agreement 
with experiment was found and this was one of the rea- 
sons to extrapolate this type of calculations to stellar 
environments, as well as to other WP nuclei. 

Measurements of the decay properties (mainly half- 
lives) of nuclei in this mass region have been reported 
in the last years [69l-l80j . The calculation of the half- lives 
in Eq. ((25|) involves the knowledge of the GT strength 
distribution and of the Qec values. In this work exper- 
imental values for Qec are used. They are taken from 
Ref. 34] or from the Jyvaskyla mass database [tI, [81| , 



when available. In Fig. |4] the measured half-lives are 
compared to the QRPA results obtained from the equi- 
librium deformations of the various isotopes. In general 
good agreement for the N = Z WP is obtained. Also for 
the more stable N = Z + 2 the agreement is very reason- 
able, except for the heavier Cd an Sn isotopes, where the 
half-lives are overestimated. The half-lives of the more 
exotic isotopes are fairly well described by QRPA. 



C. Stellar weak decay rates 

Figures [5]fTBl show the decay rates as a function of the 
temperature T. On the left-hand side (a) one can see 
the decomposition of the total rates into their contribu- 
tions from the decay of the the ground state Ogg 1+ 
and from the decay of the excited state 2+ 1+, 2+, 3+ 
in the parent nucleus. The middle panel (b) contains 
the decomposition of the rates into their /3+ and cEC 
components evaluated at various densities (pYe)- On the 
right-hand side (c) the total rates for various densities are 
presented. The gray area is the relevant range T = 1 — 3 
GK for the rp process. Each figure contains the results 
for three isotopes. The results corresponding to the more 
exotic ones are displayed on top, whereas the results cor- 
responding to the more stable isotopes appear on the 
bottom. In the middle we find the intermediate isotopes, 
which in most cases correspond to the WP nuclei. 

The results decomposed into their contributions from 
various parent states (a) show that the decay from the 
ground state is always dominant at the temperatures 
within the gray area of interest. The contributions of 
the decays from excited states increase with T, as they 
become more and more thermally populated, but in gen- 
eral they do not represent significant contributions to the 
total rates and can be neglected in most cases. Never- 
theless, there are a few cases where these contributions 
should not be ignored, which correspond to those cases 
where the excitation energy of the 2+ excited state is very 
low. This is the case of the middle-shell nuclei Kr, Sr, Zr, 
and Mo, where the contributions of the low-lying excited 
states compete with those of the ground state already at 
temperatures in the range of rp process. The effect on 
the rates of the d ecay from excited Oj states was also 
considered in Ref. [I4I in the case of Kr and Sr isotopes. 
It was concluded that in general their relative impact is 
again very small in the total rates at these temperatures. 

Concerning the competition between /S"*" and cEC 
rates (b) one should distinguish between different iso- 
topes. Thus, the more exotic isotopes appearing on the 
top of the figures show a clear dominance of the /3+ 
rates over the cEC ones that can be neglected except 
at very high densities beyond rp-process conditions. On 
the other hand, the opposite is true with respect to the 
more stable isotopes on the bottom, where the /S"*" rates 
are completely negligible. The origin of these features 
can be understood from the behavior of the phase space 
factors as a function of the available energy Qif. As it 
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was mentioned above and discussed in Ref. [17[, more 
exotic nuclei with larger Qif values favor /3+ because of 
the larger phase space factors, while the opposite is true 
for more stable nuclei with smaller Qif values. 

The interesting cases occur in the middle panels that 
correspond in most cases to the N = Z WP nuclei. Here, 
there is a competition between /3+ and cEC rates that 
depends on the nucleus, on the temperature, and on the 
density pY^. One can see that for large enough densities, 
cEC becomes dominant at any T. For low densities, 
rates dominate at low T, while cEC dominates at higher 
T, but in general there is a competition that must be 
analyzed case by case. 

Finally, the total rates in (c) are a consequence of 
the competition between /J"*" and cEC rates mentioned 
above. Since the /?+ decay rate is independent of the 
density and depends on T only through the contribu- 
tions from excited parent states, the total rates are prac- 
tically constant for the more exotic isotopes in the upper 
figures, only modulated by the small contribution from 
cEC. In the central isotopes the rates are the result of 
the competition discussed in (b), and finally in the heav- 
ier isotopes (lower figures) we can see that the total rates 
are practically due to cEC with little contribution from 
/3+. Tables containing /3+, cEC, and total decay rates 
for all the isotopes considered in this work are available 
in Ref. (sl. 



IV. SUMMARY AND CONCLUSIONS 

In summary, the weak decay rates of waiting point 
and neighbor nuclei from Ni up to Sn have been in- 
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FIG. 1: (Color online) Potential energy curves for the even-even isotopes considered in this work obtained from constrained 
HF+BCS calculations with the Skyrme force SLy4. 
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FIG. 2: (Color online) Calculated GT strength distributions for Ni, Zn, Ge, Se, Kr, and Sr isotopes obtained from their ground 
states, as well as from the shape coexisting states. The individual strengths correspond to the ground states, whereas folded 
distributions are shown for the various configurations considered in each isotope. Qec values and proton separation energies 
Sp are shown by vertical lines. 
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FIG. 3: (Color online) Same as in Fig. [2] but for Zr, Mo, Ru, Pd, Cd, and Sn isotopes. 
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FIG. 4: (Color online) Calculated QRPA half-lives compared to experimental values. 
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FIG. 5: (Color online) Decay rates (s"^)of 50.52, 54 jgotopes as a function of the temperature T (GK). (a) Decomposition of 
the total rates into their contributions from the decays of the ground and excited 2^ states, (b) Decomposition of the rates 
into their /J"*" and cEC components evaluated at different densities, (c) Total rates at various densities. The label p stands for 
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FIG. 6: (Color online) Same as in Fig. \5\ but for 56.58,60 isotopes. 
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FIG. 7: (Color online) Same as in Fig. \5\ but for 62,64,66 jgotopes. 
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FIG. 8: (Color online) Same as in Fig. \5\ but for 66.68,70 gg isotopes. 
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FIG. 10: (Color online) Same as in Fig. [5] but for '^^-''''.TSgj. jgotopes. 
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FIG. 11: (Color online) Same as in Fig. [5] but for 80'82,842j. isotopes. 
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FIG. 12: (Color online) Same as in Fig. [5] but for *'''*^'**Mo isotopes. 
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FIG. 13: (Color online) Same as in Fig. [5l but for 88.90.92j^u isotopes. 
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FIG. 14: (Color online) Same as in Fig. [5l but for 92,94,96pj isotopes. 
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FIG. 15: (Color online) Same as in Fig. \5\ but for ss.ss-ioocd isotopes. 
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FIG. 16: (Color online) Same as in Fig. \5\ but for 100.102.1043^^ isotopes. 



